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DETECTION OF ABSORPTION FEATURES IN THE X-RAY SPECTRUM OF THE NARROW-LINE QUASAR 
PG 1404+226: POSSIBLE EVIDENCE FOR ACCRETION DISK WINDS 

Surajit Dasgupta \ A. R. Rao \ G. C. Dewangan 2 , V. K. Agrawal ;i 
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ABSTRACT 

We present the results of an analysis of data from XMM-Newton and CHANDRA observations of the 
high luminosity narrow-line quasar PG 1404+226. We confirm a strong soft X-ray excess in the X-ray 
spectrum and we find rapid variability (a factor of two in about 5000 s) . When the X-ray spectrum is 
fit with a two component model which includes a power-law and a blackbody component, we find that 
low energy absorption lines are required to fit the data. If we interpret these lines as due to highly 
ionized species of heavy elements in an outflowing accretion disk wind, an outflow velocity of ~ 26000 
km s _1 could be derived. One interesting feature of the present observation is the possible detection 
of variability in the absorption features: the absorption lines are visible only when the source is bright. 
From the upper limits of the equivalent widths (EW) of the absorption lines during the low flux states 
and also from the model independent pulse height ratios, we argue that the strength of absorption is 
lower during the low flux states. This constraints the physical size of the absorbing medium within 100 
Schwartzschild radius (R g ) of the putative supermassive black hole. We also find a marginal evidence 
for a correlation between the strength of the absorption line and the X-ray luminosity. 
Subject headings: galaxies: active — quasars: individual(PG 1404+226) — X-rays: galaxies 



1. INTRODUCTION 

Narrow-line Seyfert 1 galaxies (NLSls) have very re- 
markable X-ray properties: they show evidence for strong 
excess of soft X-rays (dominant below ~2 keV) above 
the hard X- ray continuum ext rapolation and rapid X-ray 
variability l(Boller et al.lll996lK Complex absorption fea- 
tures are also c ommon i n many of these sources. Recently 
iPounds et alJ l)2003albfl detected several absorption lines 
in the narrow emission line quasars PG 1211+143 and 
PG 0844+349 which are blue-shifted indicating relativis- 
tic outflow of highly ionized material. They suggest that 
these outflows form a significant component in the mass 
and energy bud gets of systems ac creting at or above the 
Eddington rate flKing fc Poundsll200l . 

PG 1404+226 (V=15, M=-23.4, z=0.098) is one of 
the most extreme narrow line Fe II quasar. Its Hp 
FWHM is 880 km s" 1 l|Boroson & Greed Il99l . The 
R OS A T spectrum (0.1 - 2 keV) of PG 1404+226 is steep 
(r ~ 3) and shows rapid (factor of 2 in 10 hours) flux 
variability. Moreover flux selected spectral analysis re- 
vealed the presence of an absorption edge around 0.8 - 1 
keV whos e energy shift to h igher value when the source 
brightens ljUlrich fc Molendilll996j) . The X -ray spectrum 
of P G 1404+226 was variable during an ASCA observa- 
tion ijComastri et alJll997t lUlrich et al.lll99^ hereafter 
U99) and characterized by a strong soft excess below 2 
keV whose luminosity in 0.4 - 2 keV band (7 x 10 43 ergs 
s -1 ) was a factor o f 3 gre ater than 2-10 keV luminos- 
ity. iComastri et alJ l)1997j) have also found an absorption 
edge at 1.07 keV and they su ggested an over abundance 
of iron. iLeighlv et alJ (|1997fl compared the absorption 
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feature found in ASCA spectrum with absorption by ion- 
ized oxygen and derived a high velocity of ionized outflow 
from that object. 

In this Letter we present the XMM-Newton and 
CHANDRA observations of PG 1404+226. We detect 
absorption features in the XMM-Newton data and find 
some evidence for such lines in the CHANDRA data. 
We have reanalyzed ASCA data which corroborates this 
conclusion. The source shows high variability in all the 
three observations and most remarkably, the line features 
show variability in time scales of ~ 5000 s, giving a di- 
rect dynamical size constraint of 100 R g for the absorbing 
medium. 

2. DATA AND ANALYSIS 

PG 1404+226 was observed by XMM-Newton on 2001 
June 18 using the European Photon Imaging Camera 
(EPIC), and the reflection grating spectrometer (RGS) 
for about 21 ks. RGS data are not useful because of poor 
signal to noise ratio. The observation data files were pro- 
ce ssed and filtered using the same criterion as discussed 
in lDasgupta et all l)2004|) . The high energy particle back- 
ground flaring intervals were excluded. This resulted in 
'good' exposure time of ~ 14 ks for the EPIC-PN. The 
net count rate is 0.66 s _1 . The Chandra High Energy 
Transmission Grating Spectrometer (HETGS) observa- 
tion of this source was performed during 2000 July 22 
for a total duration of 80 ks. The spectrum is created 
using th order image of the source because of low counts 
collected in HEG and MEG. PG 1404+226 was observed 
by ASCA-GIS/SISon 1994 July 13-14. In this paper SIS 
data are reanalyzed. The SIS was operating in 1-CCD 
mode and the data were collected in Faint mode. Stan- 
dard criterion for good time selections (U99) have been 
applied. 

EPIC-PN light curves of bin size 500 s (background 
subtracted) of PG 1404+226 in the energy range 0.3 
- 1 keV (soft), 1 - 10 keV (hard) are plotted in Fig- 
ureQ The energy bands were chosen to separate approx- 
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Fig. 1 — (a) EPIC-PN 500s bin (left), (b) CHANDRA-ACIS 2000s bin (middle), and (c) ASCA-SIS (SISO AND SIS1 added) 5000s bin 
light curves and the corresponding hardness ratio (see text). 



imately the two spectral components - a power-law and 
a soft excess c omponent generally observed from NLSls 
<ILeighrvl H999). It is evident that the X-ray emission 
from PG 1404+226 varied strongly during the XMM- 
Newton observation. The average count rates in the low 
and high flux states (demarcated by dotted lines in Fig- 
ure [TJ are 0.26 ± 0.01 and 0.69 ± 0.02 s _1 respectively 
in soft band, 0.029 ± 0.007 and 0.042 ± 0.007 s" 1 respec- 
tively in hard. The hardness ratio (defined as the ratio 
between flux of hard band to that of soft band) is de- 
creasing with time, implying that the spectrum softens 
in the high state. ACIS and SIS light-curves (Figure Q 
show similar timescale of variability. 

Photon energy spectra of PG 1404+226 and associated 
background spectra were accumulated from the EPIC- 
PN, ASCA-SIS and CHANDRA-ACIS data. The pulse 
invariant channels were grouped such that each bin con- 
tains at least 20 counts. Data above 8 keV are not used in 
spectral analysis because background counts are compa- 
rable with source counts in this regime. All the spectral 
fits were performed with the XSPEC 11.2.0 and using 
the x 2 statistics. The quoted errors on the best-fit model 
parameters are at the 90% confidence level (Ax 2 = 2.7). 
Luminosities are derived assuming isotropic emission. A 
value for the Hubble constant of H = 50 km s _1 Mpc -1 
and a standard cosmology with q = has been adopted. 

We fitted a redshift ed powerlaw mod el with a Galactic 
N H of 2 x 10 20 cm" 2 ijElvis et al.lll989|) to the EPIC-PN 
spectrum in the energy range of 2 - 8 keV. This provides 
an acceptable fit (x /dof ~ 26/24 and T -1.72). There 
is a excess emission above 6.5 keV, but addition of a nar- 
row redshifted Gaussian line near 6.4 keV does not im- 
prove the fit (Ax 2 < 4) . Extrapolation of the best fitting 
2-8 keV power law to 0.3 keV shows a huge soft excess 
in the spectrum. Adding a redshifted blackbody gives a 
reasonable fit {x 2 /dof ~ 189/162, T ~ 1.46, kT bb ~ 108 
eV). The model (model a in Figure OJ leaves significant 
residuals in the 0.8 - 1.2 keV (observed frame) band and 
at 3 keV. Addition of Gaussian absorption lines at 1, 1.2, 
and 3 keV to the previous model (hereafter model e) gives 
a good fit with x 2 (dof) ~ 134(156) (Figure^. All the 
derived parameter values along with the errors are given 
in the Tabled We have fitted several other models to the 
data based on the previous results reported on the same 
source (Figure 0). An absorption component (model 
absori in XSPEC) along with the blackbody and power- 
law does improve the fit (x 2 /dof — 174/158 dof) but the 



absorption feature still remains (model b). Over abun- 
dance of iron (U99) cannot handle the absorption feature 
completely (x 2 /dof — 154/159, the iron abundance be- 
came more than 25 times solar abundance). Adding an 
edge with the absorption model (with super iron abun- 
dance) does not improve the fit (A^ 2 — 2). Adding an 
edge at ~1 keV with the two component model (model 
c) can handle the absorption feature at 1 keV (x 2 /dof ~ 
153/160) but the 1.2 keV feature still remains. But 
adding another edge (model d) for the — 1.2 keV fea- 
ture worsen the fit (Ax 2 — 1). One edge and one line 
also cannot handle the situation. But when we adopt 
the model with two absorption lines the fit is improved 
drastically (x 2 /dof ~ 144/158) and the residuals around 
1 keV vanished. Adding another absorption line at ~ 3 
keV improves the fit further (Ax 2 ~ 10). 

To investigate whether the flux variability is due to 
changes in the spectral parameters we have carried out 
spectral analysis at low and high flux levels. We ex- 
tracted two average spectra covering - 6.5 ks (high flux 
state) and 9.5 - 15.5 ks (low flux state) in the 0.3 - 8 keV 
band. The average EPIC-PN count rates are 0.26 ± 0.01 
s" 1 (low flux state) and 0.69 + 0.02 s" 1 (high flux state). 
We fit the data using model e. The fit parameters and 
the observed fluxes in the two flux states are given in Ta- 
ble n] The observed flux in the 0.3-8.0 keV range varied 
by a factor of 2.6 during the high and low flux states. 
The blackbody flux varied by a factor of 3 whereas the 
power-law flux varied by a factor of 1.7. The spectral 
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Fig. 2.— 0.3 - 8 keV EPIC-PN spectrum and best fitted model 
(model e, see text) and the residual spectrum (lower panel). The 
panel inside is the PHA ratio of high and low flux states. 
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Fig. 3.— 0.3 - 8 keV EPIC-PN residual spectra fitted with 
different models (see text). 



parameters are consistent with each other at 90% confi- 
dence level. Since the spectral parameters and the nor- 
malizations are strongly coupled to each other, we can 
only conclude that the soft and hard components varied 
at different ratios. The absorption lines, on the other 
hand, are not required at the low flux level. We have 
included two lines (0.99, 1.17 keV) to get an estimate of 
the upper limits to the fluxes. The 90% confidence up- 
per limits on the line fluxes during the low flux state are 
distinctly lower than those found for the high flux state. 
The upper limit on the EW has an overlap with the 90% 
confidence error for that in the high flux level. To un- 
derstand the spectral variability in more detail, we plot 
the PHA ratio of high and low flux states with energy 
(Figure [3 inner panel). While the flux difference in the 
soft emission is remarkable, the variability above 1 keV 
is much lower. There is a clear dip at around 1 keV in- 
dicating that absorption is significantly high in the high 
flux state. 

To corroborate these findings we have analysed the 
spectra from the CHANDRA and ASCA data. The de- 
tection significance above 1 keV in these two data sets is 
poor and we find a marginal evidence for the absorption 
lines in the time-averaged data. Results of the spectral 
fits to the high flux states in these two data sets (marked 
in Figure Q are presented in Table The absorption 
lines are detected in the ASCA data (Ax 2 ~ 10) and 
there is a marginal evidence for these lines in the CHAN- 
DRA data. The line energies and identific ations in case 
of AS CA data are different from that by iLeighlv et al.l 
l)1997j) probably because of different continuum model- 
ing. It is clear from Table H that the strength of the 
absoption lines increased when the flux increased. It 
is interesting to investigate whether these changes are 
correlated with each other. For this purpose we have 
divided the whole region into several parts of 1000 s du- 
ration. The energy channels were appropriately grouped 
to achieve a good signal-to-noise (a minimum of 5 counts 
per energy channel). All the spectra are fitted with the 
same model mentioned above (with C-statistics). All 
the parameters except the normalizations of blackbody, 
power-law and absorption lines are fixed to the values 
found from the fit result of the average spectrum. The 
relative normalization between the absorption lines are 
kept fixed. The variation of the free parameters with 
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Fig. 4. — Variability of normalizations (units defined in the 
footnote of Table 1) of power-law (N po ), blackbody (JV^t) and ab- 
sorption line at 1 keV (N a i, 3 ) with time. 



time is shown in Figure 0] It is clear that all the free 
parameters vary with time in a same fashion and the 
variations are similar to the count-rate variations. The 
black body flux is correlated with the line flux with a rank 
correlation coefficient of 0.79 (probability 8 x 10~ 4 ). 

3. DISCUSSION 

The spectroscopic analysis carried in this paper using 
XMM-Newton data revealed three absorption line fea- 
tures at ~ 1.0, 1.17, 3keVwithEW of 45, 60, and 364 eV, 
respectively. The line energy and EW are comparable t o 
that found in PG 1211+143 by <)Pounds et alJl2003aj) . 
We have made similar identifications for the spectral 
lines detected in PG 1404+226 and they are given in 
Table |3 For some of the lines a few possible alternate 
line identifications are also given. For the two strongest 
lines detected by the EPIC-PN (0.99, 1.17 keV) we derive 
outflow velocities of 25700±||gg and 43500±^ km s" 1 . 
However if we consider the 1.17 keV is originated from 
the line Ne IX ls-3p (1.073 keV) then the measured out- 
flow velocity will be 27000±^ km s" 1 . Though there 
is an indication of different velocities from the other line 
identifications, we cannot rule out the possibility of all 
lines originating from similar velocity structures due to 
the much lower significance level of detection of these 
lines. 

The more striking result of our analysis is the possible 
detection of line variability. From the model indepen- 
dent PHA ratio, we can see that there is strong spectral 
variability and the absorption line strength (EW) are 
also different in the two spectral states. The time scale 
of variation (~ 5000s) suggest that the absorption fea- 
tures originate in the warm absorbing material located 
within ~ 100 R g from the central source (for an esti- 
mated black hole mass of 5.2 x 10 6 M Q ). Since it is 
unlikely that the physical condition of the absorbing ma- 
terial can change with the increase in the luminosity, we 
postulate that high luminosity drives stronger winds ei- 
ther due to radiation pressure or as a result of magnetic 
reconnection. The insufficient spectral resolution of the 
instrument leads to large uncertainty of the line energy. 
Uncertainty in the central energy, shape, and possible 
multiple absorption components in a single feature may 
lead to incorrect velocities or identifications. But in spite 
of these limitations, rapid variability of line strength is 
seen. 
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TABLE 1 

Best-fit spectral parameters of PG 1404+226. 



Component 


Parameter 




EPIC-PN a 




ASCA-SIS a 


CHANDRA- ACIS a 






Average 


High 


Low 


High 


High 


Blackbody 


kT (eV) 


1143 


112+2 


119+g 




125+3 




n bb b 


9 Q7+0.09 
-0.09 


4 34+ 017 


i on + 0.90 

± -°" J -1.00 




2-62i°; 3 37 




fbb (0.3-8 kcV) e 


10.33 


14.82 


4.86 


25.82 


11.05 


Power law 


r 


1 - oy -0.12 


1 - O/2 -0.16 


^ __|_ 97 

-0.20 


. in AO 

l.3J_ 4 2 


, in A A 

l.OO_ 45 




n p l c 


4 42+ 1 ' 25 
^•^ -0.94 


5 n5 +2,30 

o.uo_ 1 70 


1 07+ 1 - 60 
± - t " -0.90 


07 1 +19.1 
cS '- 1 -14.5 


7 90+3.15 




f p l (0.3-8 kcV) c 


2.48 


2.75 


1.55 


15.32 


4.26 


Gaussian 1 


E Une (kcV) 


3 - 07 -0.16 


o n o+0.10 










n g d 


-2.80±i;|l 


-4.00±§;|f 













EW (eV) 


-348ti 9 8 * 


-438+. g!? 











Gaussian 2 


E Une (kcV) 


-, 17 +0.02 
' -0.02 


i i 7 +0.03 
J-- 1 ' -0.02 

-12.43t 384 


1.17 


1 r;Q+0.10 
1 - oy -0.08 


1 1R+ - 66 

1 - 1 °-1.16 




n 9 d 


-6.67+J;J 3 


> -4.80 


1 c 7+12.3 
10 - ' -22.2 


> -11.76 




EW (eV) 


-60+i? 

— IS 


-84 + M 


-22+? 7 


-87+11 


-17+30 

— oU 


Gaussian 3 


B H ne (kcV) 


i.oo + g-Si 


i.oo+g-21 

— u.u^ 


0.99 


1.24+°° 3 

— U.U4 


1.02+°"> 
— u . 1 u 




™ 9 d 


-13.00±2;S| 


— 93 c:t;+ 6 - 61 

ZO.dO_g 62 


> -9.30 


-35 R+ 20 ' 8 


-13 in+ 12 09 
u.iu_ 12 31 




EW (eV) 


-44l 17 


09 -16 


oq+23 
ZJ -32 


0-1+47 




Total 


/ ofcs (0.3-8 keV) e 


12.44 


16.93 


6.52 


40.25 


15.04 




f int (0.3-8 keV) e 


15.18 


20.88 


7.76 


46.58 


17.57 




L int (0.3-8 keV) f 


0.68 


0.95 


0.35 


2.11 


0.80 


x 2g 


X 2 /dof 


134/156 


101/119 


85/88 


29/38 


31/63 




Xi/dof 


144/158 


107/121 










xi 2 / d °f 


170/160 


140/123 


85/88 


33/40 


32/65 




Xiw/dof 


189/162 


162/125 


85/88 


39/42 


34/67 



Parameter values for EPIC PN, ASCA-SIS and CHANDRA-ACIS data. 

b Blackbody normalization in units of 10 — 5 X 10 39 erg s _1 /(d/10 kpc) 2 , where d is the distance. 
c Power-law normalization in units of 10 — 5 photons cm -2 s — 1 kcV - 1 at 1 keV. 
d Gaussian normalization in unit of 10 _6 photons cm -2 s — 1 
°Flux in the unit of 10 — 13 erg cm~ 2 s — 1 
'Source luminosity in the unit of 10 44 erg s _1 

g X 2 : excluding Gaussian 1, x?2 : excluding Gaussian 1 & 2, Xi23 : excluding all Gaussian lines 



TABLE 2 

Absorption lines identified in the parametric fit to the spectrum of PG 1404+226. 



EW A x 2 /dof 
(eV) 



E a , 
(kcV) 



1 00+ Tmr 

i - uu -0.02 
1 17+O.02 
r*r 1 _n no 



-46 T 9 7 18/119 

,+26 
-16 



Instrument Ei ab 
(kcV) 



Line 



velocity 
(km s ) 



9 c »7nn+ 33uu 
43 ,00+ 5900 
27000+j£gg 



1.24 
1.59 



+0.03 
-0.04 



+0.10 
0.08 



o n7 +0.06 
■ iMI -0.16 



-61" 1 



-81 
-87 



+47 
-68 

+87 
-68 

.+208 



17/117 



6/40 
5/38 



364+^1 10/154 



EPIC-PN 0.921 
EPIC-PN 1.022 
1.073 

ASCA-SIS0 1.022 
1.211 

ASCA-SIS0 1.47 
1.211 

EPIC-PN 2.62 



Ne IX ls-2p 
Ne X Lya 

Ne IX ls-3p 
Ne X Lyct 
Ne X Ly/3 

Mg XII La 
Ne X Ly/3 

S XVI Lya 



+8800 
-11800 
7400 
7200 
,+20000 
16000 
94000+ 20500 

51 500+ 6900 
olouu - 18300 



64000] 
7200^ 
24500^ 
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